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On a Singular Integra differential Equation 
arising from a Linearised Free Surface Problem 

Patrick McLean* 
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Abstract 

A problem of linear surface waves discussed by Forbes [|Tj initially 
gave rise to a singular integrodifferential equation over the real line. 
We have been able to transform this integrodifferential equation into 
a linear second order differential equation whose solution has been 
found explicitly in terms of the sine and cosine integral functions. 
Furthermore, we have been able to recover known results of physical 
interest using more simple techniques. 
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1 Introduction 

In this paper we consider the singular integrodifferential equation 

u(x) + F 2 Hu'(x ) = f(x), —oo < x < oo, (1) 

where H is the Hilbert transform, F > 0 is a constant and f(x) is a certain 
algebraic function. We show that m is equivalent to a linear second order 
differential equation which can be explicitly solved in terms of the sine and 
cosine integral functions. 

This equation arises from a free surface problem in fluid mechanics and 
this background is discussed in section 2. In section 3 we use certain prop¬ 
erties of the Hilbert transform to convert CD to a differential equation the 
solution of which is then evaluated in terms of the sine and cosine integral 
functions. We are also able to extract certain properties of the surface. Fi¬ 
nally, in section 4 we present graphs of the solution u(x) for a variety of 
F. 


2 Background 


In |5 3 and (6. Forbes considers surface waves caused by flow over a (line) vor¬ 
tex. He models the situation using a complex potential $(z) beneath a free 
surface y = S(x). This results in a nonlinear free surface problem. Lineari¬ 
sation (assuming the vortex strength is small) and employing the boundary 
integral method results in the following singular integrodifferential equation 


(tyi 

dx 


(x,Q) + 



d 2 4*i 

d 2 t 


M) 


dt 

t — x 


1 

7r(x 2 + 1) ’ 


—oo < x < oo, 


( 2 ) 


where 0i is the real part of the linearised complex potential and the integral is 
to be interpreted in the Cauchy principal value sense. This equation appears 
as equation (2.8) in (Hj. 

If we define the Hilbert transform of a function f(x) as in [3} §15] 

urn = - /” F«L dyt (3) 

7T 7-00 V - X 


then we can express equation © in the form © where 

«W = ^fe°), f ( x) = - { -L- l (4) 
The linearised surface and complex potential can be recovered from u(x). 
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Equation 0 is an integrodifferential equation, thus to solve it uniquely 
we need to specify an extra condition. Namely, we seek a solution to equa¬ 
tion 0 which vanishes towards — 00 , that is, we shall impose the boundary 
condition 

lim u(x) = 0. (5) 

X —► — OO 

3 Analytical Solution 

In this section we convert the integrodifferential equation 0 to a differential 
equation from which we obtain various results regarding the solution u(x). 
In particular, we obtain a representation of the function u in terms of sine 
and cosine integral functions, we determine the asymptotic behaviour of u(x) 
towards +00 and we evaluate w(0) and w'(0). 

3.1 Conversion to a differential equation 

We can convert the integrodifferential equation 0 into a second order dif¬ 
ferential equation. The resulting equation can be solved and the boundary 
condition 0 invoked to give a representation of u(x) as an integral. 
Differentiating 0. applying H and multiplying by F 2 gives 

F 2 Hu' + F a H{Hu'Y = F 2 Hf. (6) 

Using the results H 2 u = —u, see 0 §15.1 (2)] and Hu' = (Hu)', see (JJ §15.1 
(8)] we have that H(Hu')' = —u". Thus, 0 becomes 

F 2 Hu' - F 4 u" = F 2 Hf. (7) 

Subtracting © from 0 gives the 2nd order differential equation 

FV + u = f- F 2 Hf. (8) 

The general solution to the homogeneous part of this equation is 

Asin^+Fcos^) (9) 

for arbitrary constants A and B. By the method of variation of parameters, 
see, for example, [3j, we find that a particular solution, u p say, which is zero 
at — 00 is given by 

M x ) = J sin (-pT^) [/(*) - F2 Hf(t )] dt. (10) 
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Thus, the general solution of equation (jHJ) is given by 

u(x) = Asin(^) + B cos(J^) + u p {x). (11) 

Imposing condition © we have, since u p (— oo) = 0, that both A and B must 
be zero. So, the required solution is 

u(x) = -^2 J sin(^^) [f(t) - F 2 Hf(t )] dt. (12) 

Since Hf = ( Hf)' we can apply integration by parts to give 


sin f)'(t)dt = 


F 2 






cos( — )Hf(t)dt. 


F 2 


The following Hilbert transform is to be found in |1J §15.2 (10)] 




-) = 


-x 


x 2 + 1 x 2 T 1 ’ 


from which we see that 


lim Hf(t) = 0. 

t —► — OO 


From m, m and (USD we have the representation 
IF 


u(x) = 


F 2 


x /. . . .x t . , . 

sm( -^r- )f(t ) + cos (-^ r )Hf(t) 


F 2 


dt. 


(13) 

(14) 

(15) 

(16) 


3.2 Trigonometric Integrals 


We follow the standard reference work E §5.2] and define the trigonometric 
integrals as follows. For z Ed! define the sine integral Si( 2 ) as 


Si(*) = 


sin t 
t 


dt. 


(17) 


The sine integral Si is an entire function. A commonly used notation is 


si ( 2 ) = f — Si( 2 ). 

For z E(P such that | arg^| < 7 r define the cosine integral Ci( 2 ) as 

rz cos t — 1 


Ci( 2 ) = 7 + log 2 + 


'0 


t 


-dt. 


(18) 


The cosine integral Ci has a branch cut discontinuity along the negative real 
axis. 
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The functions Si and Ci occur as Fourier sine and cosine transforms. 
Specifically, for | arga| < n and y > 0 we have that 



cos (xy) 
a + x 
sin (xy) 
a + x 


dx 

dx 


— si (ay) sin(ay) — Ci(ay) cos (ay) 
Ci (ay) sin(ay) — si(ay) cos (ay), 


see m § 1.1 ( 9 )] and 0 §1-2 (10)], respectively. 


(19) 

( 20 ) 


3.3 Evaluation of Solution in Terms of Trigonometric 
Integrals 

In this section we evaluate u(x) in terms of the trigonometric integrals Si 
and Ci. 

Invoking the change of variable x — t — s in equation m we obtain 


u(x) 



sin (™) 


F 2 (x 


1 

s) 2 + 1 


+ 


x — s 
- s) 2 + 1 


ds. 


( 21 ) 


On making partial fraction expansions of the algebraic expressions in 
equation m and using the identities PI). ©, together with standard 
trigonometric identities we obtain the representation 


Sill 


ulx = 


F 2 ■ 


7 rF 2 e^ 

| cos(^) 

7 iF 2 e'p I 


1 

2 i 


Ci(— 

Ci(- 


F 2 


x 


Ci(- 


F 2 


+ Ci(- 


x + i 

p2 ) 

x + i 


F 2 


1 

~ 2 

1 

+ 27 


QV X ~F , Q .f x + ^ 

Si(—™~) + Si(- 


Si(- 


p2 > v p2 

x-i. n .. x + i. 

)- Sl (- ft) 


+ 


F 2 


F 2 


This representation allows evaluation of u(x) for arbitrary x in terms of 
the functions Si and Ci at complex arguments. High precision evaluation 
of these functions is provided by computer packages such as Mathematica. 
For Fortran code for evaluating the sine and cosine integrals at complex 
arguments see [2j. 

An alternative method of evaluating u(x) is to view m as a linear 
combination of the Fourier sine and cosine transforms 



s) 2 + 1 


ds 



cos (t^) 


F 2 '(x 


x — s 
- s) 2 + 1 


ds , 


(23) 


and use the method of jjj to numerically evaluate them. 
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3.4 Asymptotic Behaviour of Solution 


In this section we examine the asymptotic behaviour at +oo of the solution 
m- Rewrite m as 


u{x) = 


F 2 J 
1 


x — t \ „ . \ ,x — t \ „ . \ 

- cos (——)Hf(t) 


—oo _ 

PCO 


p2 W v 

- cos 


dt 


dt. 


(24) 

(25) 


The second integral in this representation vanishes as x —> oo, so we have 
that 


u(x) ~ 


1 

F 2 


p0 ° . ,x-t. r , N , 1 

sin( — )f(t)dt — 


F 2 


F 2 


cos( X „„ )Hf(t)dt , (26) 


F 2 


for large positive x. We shall need the following two properties of the Hilbert 
transform: 



y)Hg(y)dy = 



y)g(y)dy. 


(27) 


and 

Ft cos(x) = — siu(x), (28) 

which are to be found in [3 §4.2 (16)] and |U §15.2 (47)], respectively. 

Thus, (BED becomes 

1 T — 1 1 /*°° 7 * — / 

«(t) ~ ^7 J J Hcos(-y r ) f(t)dt (29) 

9 r°° r — t 

= J ~ 2 J jin(—)f(t)dt. (30) 

Now, using the fact that the sine function is odd and / is even it follows that 


u(x) 


nF 2 


sin (^)/ oo COS (^ 7 ) 


t 2 + 1 


dt. 


(31) 


On using the Fourier cosine transform to be found in [41 §1.2 (11)] 



cos (xy) 
x 2 T a 2 


dt = 


7T 

2a 


,-ay 


a,y> 0 


(32) 


we obtain 


u(x) ~ 2e" 


i 

'jv 


Sill 


X , 

~p2‘ 


X 


oo. 


(33) 


This agrees with [3 eqn 3.4], Thus, we have found the asymptotic behaviour 
at Too of the function u(x). That is, the solution oscillates sinusoidally at 
Too. 
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3.5 Behaviour at origin 

We now note the value of u(x) and its derivative at the origin x — 0. Evalu¬ 
ating m at x = 0 we have 


u(0) 


1 r- sim'jJ,) - s':ns(fi) 

- —1 - “ S 


(34) 


The following identities appear in [5j as equations 2.2 (14) and 1.2 (12) 
respectively 


f°° sin (xy) 

Jo x 2 + « 2 

f 00 x cos (xy) 
J o x 2 + a 2 


dx 

dx 


[e -ay Ei(ay) — e ay Ei(—ay)] 

-T [e-““ Ei(ay) + e ay Ei(-oj)] , 


(35) 

(36) 


where a, y > 0 and the exponential integral function Ei is defined as 



(37) 


see [I,; eqn 5.1.2]. Thus, we are able to express w(0) as 


e "f 5 ' 1 

“(°)=^FEi (j5 ). 


Similarly, we can derive that 


u'{ 0) = 


_ i 

e f 7 

~T r 


(38) 


(39) 


4 Linearised Surface Profile 


In this section we present graphical representations of the surface prohle for 
three Froude numbers F. From equations (2.1) and (2.5a) of |Sj the surface 
prohle is given by 

S(x) = -eF 2 ^jjj-(x, 0) + 0(e 2 ). (40) 

We take e = 1, F = 0.1,1,10 and plot u(x) in figure |U We note that 
in each case u(x) is a superposition of an impulse due to the vortex and a 
sinusoidal wave downstream. Qualitatively the main differences are scale and 
the relative strengths of these two contributions. 
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Figure 1: Surface profiles for F = 0.1,1,10 
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